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ABSTRACT 



We present the analysis of the broadened, flourescent iron Ka line in simultaneous XMM-Newton and RXTE data from the black 
hole Cygnus X- 1 . The XMM-Newton data were taken in a modified version of the timing mode of the EPIC-pn camera. In this mode 
the lower energy threshold of the instrument is increased to 2.8 keV to avoid telemetry drop outs due to the brightness of the source, 
while at the same time preserving the signal-to-noise ratio in the Fe Ka band. We find that the best-fit spectrum consists of the sum 
of an exponentially cut off power-law and relativistically smeared, ionized reflection. The shape of the broadened Fe Ka feature is 
due to strong Compton broadening combined with relativistic broadening. Assuming a standard, thin accretion disk, the black hole is 
close to rotating maximally. 
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1. Introduction 

Relativistically broadened lines from accreting black holes have 
an important diagnostic potential because they allow us to mea- 
sure directly properties of these black holes such as their spin, 
the inclination of the accr etion disk, and the emissi vity profile 
of the accretion disk (e.g.. lRevnolds & Now ak 2003). While by 
now a large number of high signal-to-noise ratio (S /N) mea- 
surements of the line shape and p arameters exist for many 
Galactic blac k holes (iMillerl I20071) and for Activ e Galactic 
Nuclei (e.g.. lMiniutti et al.ll2007t iFabian et al.ll2009l) . high S/N 
measurements of the line profile and black hole parameters of 
the canonical black hole, Cygnus X-1, are still lacking. This 
lack is mainly due to a combination of source brightness, which 
causes the measurement to be domi nated by systema tic effects 
in the detector calibration (see, e.g.. lMiller et al.ll2010t) . as well 
as a complex X-ray spectrum, which necessitates broad band 
(3-1 00 keV) mea s ureme nts from multiple satellites. See, how- 
ever, iMiller et al.1 (|2002) for recent work on the Fe line region 
with Chandra's gratings, and lNowak et alj (1201 ll and references 
therein) on joint Suzaku/RXTE measurements of the broad band 
spectrum. 

In this letter we present an analysis of the broad band spec- 
trum of Cygnus X-1 using data from high S/N observations 
of the Fe line band from XMM-Newton and simultaneous, 3— 
120keV RXTE data. In Sect. [2] we present the modified timing 
mode, a special form of the XMM-Newton timing mode, ded- 



icated to bright sources with fluxes up to a few lOOmCrab. In 
Sect.[3]we discuss the data reduction in greater detail. Section|4] 
presents the analysis of the relativistically broadened line and 
the reflection continuum. 

2. The EPIC-pn modified timing mode 

To determine the best Fe line parameters possible, very high 
S IN observations are required in a time that is as short as pos- 
sible to avoid complications due to source spectral variability. 
The best instrument currently fulfilling thes e requirements is th e 
EPIC-pn camera onboard XMM-Newton (Stru der et at]l200ll) . 
This camera has two modes that provide the fast rea d-out capa- 
bility required to avoid pile-up in bright source data (Ness et al. 
2010). Here, the term "pile-up" summarizes both "energy pile- 
up", where more than one X-ray photon impacts the same CCD 
pixel during one CCD read-out cycle, and "pattern pile-up" or 
"grade migration", where multiple X-ray photons hit adjacent 
CCD pixels during one read-out cycle and the resulting charge 
distribution in the CCD is wrongly interpreted as originating 
from a single, higher energy, X-ray. For satellites with moderate 
spatial resolution such as XMM-Newton, the latter effect is more 
pronounced. In EPIC-pn's burst mode, fast exposures and shifts 
are performed, interrupted by longer periods for data read out. 
The livetime of this mode is only 3%. In the timing mode, con- 
tinuous read out is available with a 99.5% livetime. This mode 
is well suited for observations with a low pile-up fraction for 
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Fig. 1. XMM-Newton and RXTE lightcurve of Cygnus X-l on 
2004 November 20/21. Here, and throughout the remainder of 
the letter, EPIC-pn data are shown in blue, PCA (PCU2) data in 
red, and HEXTE in orange, with resolutions of 100 s, 96 s, and 
40 s, respectively. 
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Fig. 2. a Fit of a simple power law to the 4.0-9.5 keV data, ig- 
noring the 5.0-8.0keV band (light gray area), b A broad excess 
due to the broad line, is visible. Note the energy shift between 
the PCA and the EPIC-pn. c Residuals after applying the gain 
shift correction. 



fluxes up to ~150mCrab, although the EPIC-pn telemetry allo- 
cation limits this mode to observations of sources of fluxes be- 
low -lOOmCrab. With ~300mCrab, Cyg X-l can therefore not 
be observed with the standard timing mode. 

In order to address the telemetry li mitation we have propose d 
a modification of the timing mode (Kendzio rra et alj 12004). 
In this modified timing mode, the lower energy threshold for 
telemetered events is increased from 0.15 keV to 2.8 keV, such 
that all events in the Fe line band can be transferred without 
telemetry drop outs, and the telemetry allocation to the EPIC-pn 
is maximized by switching off the EPIC-MOS cameras, which 
are piled-up anyway for sources above ~35 mCrab. 

EPIC data from events in which the charge cloud produced 
by an X-ray is distributed over more than one CCD pixel, 
so-called "split events", are only recombined on ground. In 
the modified timing mode, split partners with energies below 
2.8 keV are not telemetered, which results in a slight degrada- 
tion of the energy resolution and necessitate s the generation of a 
dedic ated respo nse matrix. As discussed by Wil ms et alj ((2006a) 
and lFritzl <H.K>8). this effect can be fully taken into account by 
measuring the energy dependent probability distribution of split 
events from archival Timing Mode observations and modify- 
ing the official Timing Mode response matrix accordingly. A 
full de scription of the ca libration will be given in a forthcoming 
paper dDuro et alj|201 ll) . The XMM-Newton data analysis and 
response matrix generation was based on the Science Analysis 
Software (SAS), version 10.0.0, and the newest calibration files. 

It has recently been found that so me Timing Mode observa - 
tions are affected by "X-ray loading" dDone & Diaz Tr igo 2009), 
i.e., the contamination of EPIC-pn's offset map by source X- 
rays. Since the offset map is subtracted from the EPIC data be- 
fore they are telemetered to ground, loading distorts the X-ray 
spectrum and effectively introduces a CCD-column dependent 
lower energy threshold. In the modified timing mode, the offset 
map is determined with the EPIC-pn filter in the blocking po- 
sition. As no source X-rays enter the offset map, the mode is 
unaffected by X-ray loading. 



3. Data reduction 

3.1. Data reduction 

In 2004 November and December XMM-Newton performed 
four modified timing mode observations of Cyg X-l. Here, we 
concentrate on the simultaneous XMM-Newton and RXTE ob- 
servations performed 2004 November 20 and 2 1 {XMM-Newton 
ID 0202760301), during which Cyg X-l was in a transitional 
state between th e hard state and the thermally dominated soft 
state (Fritz 2008). We concentrate on this 17.4ks long observa- 
tion since it shows the least variability on timescales of minutes 
to hours (Fig.Q]) and thus provides the cleanest X-ray spectrum. 
Result s from all observ ations will be presented in a subsequent 
paper dDuro et al.ll201 lb . The increased soft X-ray emission dur- 
ing the intermediate state potentially leads to pile-up in the cen- 
ter of the point spread function. For the present study, we there- 
fore ignore the innermost three columns. In order to avoid any 
transient effects due to the 2.8 keV lower energy threshold, we 
limit ourselves to the 4.0-9.5 keV band. 

The simultaneous 6.18 ks of RXTE data (observation IDs 
90104-01-02-00, -01, and -03) were reduced with HEASOFT ver- 
sion 6.9. We used th e 4-40 keV data f rom the Proportional 
Counter Array (PCA; Uahoda et al.ll2006l) and the 20-112keV 
data from the High En ergy X-ray Timing Experiment (HEXTE; 
Rothschil d et al.lTl998l) and ignored data taken within 10 min- 
utes of passages through the South Atlantic Anomaly and dur- 
ing times of high particle background. All spectral analysis 
was done with the Interact ive Spectral Interpretati on System 
dHouck & Denicolal l2000l iNoble & Nowakl l2008h . HEXTE 
spectra were rebinned to S /N = 30, EPIC-pn to S/N = 100, 
which resulted in a resolution of 40 eV at 6.4 keV and oversam- 
pled the detector response by a factor ~4. 

3.2. Simultaneous XMM-Newton and RXTE Fits 

The 3-120keV spectrum of Cyg X-l can be well described as 
the sum of a soft excess with a temperature of <500 eV, an ex- 
ponentially cutoff power law with a cutoff at a few lOOkeV, and 
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a Compton reflection component dNowak et alJ 1201 ll and ref- 
erences therein). Alternatively, a variety of other, more physics 
based, components h ave been used to describe the continuum 
dNowak et alJ 1201 I t iMalzac et al.1 120081: IWilms et al.1 I2006H: 
Markoff et al. 2005), but these are virtually indistinguishable 
from this empirical model. In addition to the continuum com- 
ponents, Chandra -HETGS observations show the presence of a 
narrow line at 6.4 keV emitted by neutral i ron in material sur - 
rounding the black hole at farther distances dHanke et al.l l2009). 

As shown in Fig. [2] a simple continuum model is not suf- 
ficient to describe the data. Strong positive residuals in the re- 
gion between 5 and 7 ke V remain, which are due to the broad Fe 
Ka line. However, a slight mismatch in the RXTE and XMM- 
Newton residuals is also apparent (Figure [2J}), with the EPIC 
line appearing to be shifted to a slightly higher energy. This shift 
is confirmed by measuring the energy of the 6.4 keV narrow line 
feature, which is found at 6.6 keV, inconsistent with all available 
Chandra measurements for this feature in Cyg X- 1 . 

We attribute this difference to the treatment of charge trans- 
fer efficiency (CTE) effects in the SAS. When observing bright 
sources with a CCD such as the EPIC-pn, traps in the sili- 
con crystal can be saturated with source electrons, which im- 
proves the CTE. While the latest SAS versions take this effect 
into account, their calibration breaks down for very high count 
rates. An overcorrection of CTE effects then leads to inferred 
photon energies that are too high. We take this effect into ac- 
count by correcting the assigned photon energies with a model 
of the form E K!i \ = £obs/s + A£. The model parameters s and 
AE are determined using measured energies of narrow spectral 
components or from continuum fi tting. As will be discussed in 
greater detail bv lDuro et al.1 (1201 ll) . Chandra's HETGS observa- 
tion (obsid 3814) was taken at the same flux level as our XMM- 
Newton/RXTE observations. In addition to the 6.4 keV line, 
Chandra reveals two Fe Ka absorption lines due to Fe xxv and 
Fexxvi at 6.646 keV and 6.955 keV. Fits show that full consis- 
tency between the modified timing mode spectrum and HETGS 
is obtained for s = 1 .02 and AE = ke V. This overcorrection of 
the EPIC energies by 2% is in line with typical CTE corrections 
for bright sources. An initial fit to the simultaneous 4-10keV 
EPIC and PCA data agrees with this result (Fig. |2t), showing 
that the CTE correction can also be obtained from continuum 
fitting. The remaining differences between the EPIC and PCA 
data are due to narrow features that are not resolved by the PCA. 

4. The relativistic broadened Fe Ka line of Cyg X-1 

We now turn to modeling the broadened Fe line. To describ e 
its profile we use the relline-model of Daus er et a D d2oloh . 
a relativistic line model for thin disks around black holes with 
spins -0.998 < a < +0.998, where a < indicates that the 
angular momenta of the black hole and the disk are antiparallel. 
We assume a radius-dependent line emissivity per disk unit area 
that scales oc r~ e , where e = 3 for a thin accretion disk. The disk 
extends from the (a dependent) innermost stable circular orbit to 
400GM/c 2 . 

As expected, modeling the 4-10keV spectrum with a sim- 
ple relativistic line and a power-law continuum results in a bad 
fit (xLa ~ 1 -5) an d unphysical line parameters as this approach 
does not take reflection self-consistently into account. We there- 
fore extend the data to the full 4-120keV band provided by 
EPIC-pn, PCA, and HEXTE, and model reflection of an ex- 
ponentially cutoff powe r law self-consistent ly using the reflec- 
tion model reflionx (Ross & Fabian] 120051) . This model in- 
cludes transitions from the most important ions, including Fe vi 
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Fig. 3. a Unfolded XMM-Newton and RXTE data and best- 
fit model components. Green line: reflection continuum in the 
frame of the disk (note the strong Compton broadening; the "jit- 
ter" above 50 ke V is due to the numerical resolution of the reflec- 
tion model). Purple line: relativistically smeared reflection com- 
ponent, b Measured count rate spectra, c Residuals of the best-fit 
setting the relativistic convolution to zero, d Best-fit residuals in- 
cluding relativistic convolution. 

to Fe xxvi. In order to account for the relativistic smearing we 
convolve the refl ected spectrum with a relativistic convolution 
model (relconv; Dauser et al. 2010). We assume that the angu- 
lar distribution of the fluorescent photon s is isotropic and du e 
to photons emerging from a hot corona dSvoboda etalj [2009). 
We also included narrow Ka absorption lines from Fe xxv and 
Fe xxvi and a soft excess in the model. 

We first modeled the continuum using the canonical value of 
e = 3 for a thin accretion disk. As shown in Fig. [3] and Table Q] 
this model describes the 4-120 keV energy spectrum well, with 
a reduced x^ ed - 1.10. The line profile indicates a high angu- 
lar momentum of the black hole (a = 0.88^ ^, Fig. gp), with 
additional line broadening caused by strong Compton broad- 
ening in the strongly ionized reflector (ionization parameter 
£ = 1400^jy[j erg cm s _1 ), which smears out all other discrete fea- 
tures such as edges in the reflection spectrum (Fig. [3] green line). 
The inferred inclination is 32° + 2°, which agrees well with most 
estimates for the in clination of the system (e.g ., 32° < i < 40°, 
Nink ov et al J 19871 i < 55 °. ISowers et al.lll998L or 36° < i < 67° 
Davis & Hartmann 1983. The Fe abundance is slightly higher 
than solar. 
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Table 1. Best fit results to gainshifted EPIC-pn, PCA, and 
HEXTE data assuming an exponentially cutoff power law, a 
narrow emission line at 6.4 keV, two resonant absorption lines, 
disk black body emission, and relativisticaly convolved reflec- 
tion. Fluxes and normalizations, A,, are determined relative to 
the RXTE-PCA. cepic and chexte are flux normalization con- 
stants for the EPIC and HEXTE, respectively. Uncertainties are 
at the 90% level for one interesting parameter. 
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e free 


6 frozen 


r P . 

£foid [keV] 


1 663 +tim 

j J7+0.09 

29Q+%8 8 


1.670 ±0.018 

i 00+0.07 
-0,08 


A hb 


{0M$$xVfi 


(0.54 ± 0.07) x 10 4 


^6 4kev[cgs] 


(0.6 ± 0.2) x 10- 3 


(0.69+°;^) x 10- 3 


BFexxv Ka [keV] 


6.646 


6.646 


T\ 


(0.5 ± 0.3) x 10- 2 


(0.4 ± 0.3) x 10~ 2 


BFexxvi K„ [keV] 


6.955 


6.955 
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Fig. 4. ax 2 significance contours for the emissivity index, e, and 
the spin, a. Contours are based on A^ 2 = 2.30, 4.61, and 9.21, 
i.e., 68%, 90%, and 99% confidence for two parameters of inter- 
est, with respect to the best fit case with e let free. Two minima 
are apparent, one indicating a low-spin black hole at an unphys- 
ically large e, and one indicating an almost maximally spinning 
black hole consistent with emission from a thin accretion disk 
(e = 3). bx 2 as a function of a for e — 3. 



is consisten t with spin mea surements from the accretion disk 
continuum dGou et al.ll201 lb . Being in a high mass X-ray bi- 
nary, the black hole in Cyg X-l is comparatively young. Our 
result indicates that it must have been born at close to maxi- 
mum spin as expected in typical scenarios for black hole for- 
mation in supernovae (iMacFadven & Wooslevlll999l) . In a next 
step, all four modified timing mode observations will be consid- 
ered dDuro et al.ll20Tlb . 



In order to check the uniqueness of this thin disk solution, 
we performed where e was let free. Figure |4] shows signifi- 
cance contours for this case. Most best-fit parameters are con- 
sistent between both fits, specifically, the inclination remains 
around 35° and the reflector remains strongly ionized. The emis- 
sivity, however, increases to e = 10 (the hard limit) while a 
decreases to 0. Broad lines can therefore be generated either 
by strongly concentrating all available emissivity to the inner- 
most regions of a disk around a Schwarzschild black hole, or 
by emission from a standard disk around a maximally rotating 
black hole. Both solutions are comparable in statistical qual- 
it}fl Given that most physical scenarios for very s teep emissivity 
profil es such as strongly torqued accretion disks (lAgol & K rolik 
2000) or strong illumination of the innermost disk due to strong 
light bending from a hig h latitude source ("lamppost models"; 
Martocchia & Matt 1996) also require high a, the high spin so- 
lution is preferred on physical grounds. Cygnus X-l is therefore 
plausibly close to rotating maximally. 

To conclude, this first analysis of a data set taken in the 
EPIC-pn modified timing mode confirms the capability of the 
EPIC-pn camera to observe sources up to several lOOmCrab, 
yielding high S /N measurements of the Fe line band. Coupled 
with broad band data from other satellites such as RXTE, 
INTEGRAL, or Suzaku, the modified timing mode allows the 
high precision measurement of the physical parameters of black 
hole candidates, including their inclination, spin, and emissivity 
profile. Applying this analysis to Cyg X-l shows that its black 
hole is likely a Kerr black hole if one assumes the accretion disk 
to be flat and illuminated from a Compton corona. This result 

1 For Ax 2 = 7 and 235, respectively 236 degrees of freedom, the F- 
test indicates no improvement by introducing a free disk emissivity at 
the 99% level. 
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